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Abstract 

In the template-assistance model, normal prion protein (PrPC), the pathogenic cause of prion 
diseases such as Creutzfeldt- Jakob (CJD) in human, Bovine Spongiform Encephalopathy (BSE) 
' in cow, and scrapie in sheep, converts to infectious prion (PrPSc) through an autocatalytic 

process triggered by a transient interaction between PrPC and PrPSc. Conventional studies 
suggest the S1-H1-S2 region in PrPC to be the template of S1-S2 /3-sheet in PrPSc, and the 
conformational conversion of PrPC into PrPSc may involve an unfolding of HI in PrPC and its 
. refolding into the /3-sheet in PrPSc. Here we conduct a series of simulation experiments to test 

^ ' the idea of transient interaction of the template-assistance model. We find that the integrity 

CO , of HI in PrPC is vulnerable to a transient interaction that alters the native dihedral angles at 

■ residue Asn^'*'^, which connects the SI flank to HI, but not to interactions that alter the internal 

' structure of the SI flank, nor to those that alter the relative orientation between HI and the S2 

2 ! flank. 
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Prion protein (PrP) in its infectious form is the pathogen that causes several prion diseases such as 



. Creutzfeldt-Jakob (CJD) in human. Bovine Spongiform Encephalopathy (BSE) in cow, and scrapie 

I in sheep |30| . Two reviews recently summarize past studies and discuss molecular mechanisms of 

the prion disease to understand physiological function of PrP and pathogenic pathways [H [2] . 

An in vitro experiment conducted by [7j provides strong evidence for the protein-only hypothesis. 
Fig. [1] shows an NMR structure of the C-terminal of mouse PrP in its native form (PrPC) (PDB 
code: 1AG2). It contains 103 residues from Gly^^'' to Tyr^^^ classified into secondary structures 
and surface loops [3T]. These include three a-helices: HI, residues 144 to 152, H2 (173-193), and 
H3 (200-216); two /3-strands: SI (129-131) and S2 (161-163), which form an anti-parallel /3-sheet; 
six loops: LI (124-128), L2 (132-143), L3 (153-160), L4 (164-172), L5 (194-198), and L6 (217-226). 
In what follows we shall refer to the S1-L2 segment as the SI flank, or Fl, the L3-S2 segment as 
the S2 flank (F2), and the segment from SI to S2, inclusive, as the S1-H1-S2 peptide. Homologues 
of PrP in other organisms generally have residue numbering that differ from the mouse numbering 
given in Fig. [TJ unless explicitly otherwise specified, in this text we will use the mouse numbering. Fig. 1 is here. 

*This paper has been accepted for pubcliation in European Biophysical Journal on Feb 2, 2009. 
^Present address: Department of Oncology, University of Alberta, Edmonton, AB T6G 1Z2 Canada. E-mail: 
chih-yuan.tseng@ualberta.ca; richard617@gmail.com 
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Experimental investigations suggest that the pathogeny of prion diseases is characterized by the 
unfolding of PrPC followed by misfolding into the infectious scrapie isoform (PrPSc) [22] , which 
involves conformational changes in the C-terminal residues 121-231 but no chemical reaction [roill8| . 
Muramoto et al. showed the H2 and H3 helices seem to be stabilized by disulfide bonds and likely 
to have the same conformation in PrPC and PrPSc f27]. They also found the deletion of both H2 
and H3 from PrPC does not stop its conversion to PrPSc. Eghiaian et al. showed the epitope of 
antibody to be conserved in the H2 and H3 regions, again suggesting that these two helices are 
conserved during conversion [13]. 

Because the H2-H3 region seems to be conserved during the PrPC to PrPSc conversion, recent 
experimental and computational investigations have focused on the S1-H1-S2 region (residue 124 
to 163), and these have revealed important features of the prion disease. A predominantly helical 
propensity in the HI region is demonstrated by several in silico studies under different experimental 
conditions by [33j and [31] . On the question of possible mechanisms that trigger the conformational 
conversion, in vitro studies suggest that altering the pH level of the solvent, which varies static 
electric interactions, might destabilize the HI helix and trigger conversion [39l[3. Similar effects 
were observed in in silico studies by [21] , [22j and by [9]. In simulations on mouse prion 1AG2, 
Guilbert et al. found that a major modification of dihedral angles around the residue 125 was 
required for the formation of a /3-strand in residues 121-133 [15]. In addition, DeMarco and Daggett 
show detachment of the S1-H1-S2 region from the rigid H2-H3 region to accommodate misfolding 
and the HI helix is only partially unfolded TU]. By studying the isolated S1-H1-S2 peptide, Sharman 
et al. and Ziegler et al. both showed a predominantly helical propensity in the HI region for PrP in 
water as well using CD and NMR methods [3S1 |3D] . This property is further demonstrated by Dima 
and Thirumalai's [H] and by Watzlawik et al.'s studies [37. When Watzlawik et al. experimentally 
clarified role of the HI helix in the conversion, they found the HI helix does not convert into (3- 
sheet yet it promotes aggregation. In addition, Kozin et al. found that residues 142-166 in human 
numbering (same numbering as in mouse) has the propensity to form a /3-hairpin around residues 
153-156 at neutral pH level as well [20], which they believe to be the event that drives the conversion. 
In addition, Derreumaux showed the residues 127-164 region has two equi-energetic conformations 
with /3 or a features |llj . Megy et al. provide an experimental evidence for this structural duality 
[2S]. They found the /3 structure in the region spanning residues from 142 to 167 very similar to 
NMR spectroscopy f3 structure. These results may be summarized as follows: in the PrPC to PrPSc 
conversion the S1-H1-S2 region plays an important role while the H2-H3 region plays at most a 
passive role, and during the conversion the /3-sheet breaking in PrPC is likely to be the first barrier 
[3] and detachment of the S1-H1-S2 region from the rigid H2-H3 region is another barrier [TOj . 
The HI helix is only partially unfolded. Thus, here we only aimed to study HI stability using the 
S1-H1-S2 peptide. 

One of the models advanced for the pathogenesis of the prion disease is the template-assistance 
model [30^ llT] [36] . In this model it is assumed that PrPC, normally more stable than PrPSc in 
isolation, would in the presence of PrPSc convert to the latter via a transient catalytic interaction 
with it. The implication is that a dimer of PrPSc's is energetically more stable than a system of non- 
interacting PrPC and PrPSc, which was confirmed by Morrissey and Shakhnovich's computational 
analyses [26]. When there are other PrPC present, the initial autocatalytic process would then lead 
to a propagation of PrPC to PrPSc conversion. Because PrPSc's always appear in aggregated state 
and not in isolation, its structure is not precisely known at present [9l[20], rendering an investigation 
of the conversion-causing transient interaction between PrPC and PrPSc problematic. Nevertheless, 
an in-principle feasibility of the template-assistance model was demonstrated by Malolepsza et al. 
[24j . In their computer simulations the PrPSc was approximated by peptides with /3-sheet structure, 
and the authors found that such peptides were able to induce conversion of peptides with a-helix. 

A key assumption of the template-assistance model is that the the PrPC to PrPSc conversion 
is triggered by a transient interaction, as opposed to, say, a series of slow-acting contacts. Here, 
we use computer simulation to explore possible consequences of transient interactions that may 
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trigger the PrPC to PrPSc conversion, without exphcitly including the latter in the simulation. 
In practice, we investigate what sudden changes to the conformation of PrPC would destabilize its 
native structure. If a conversion triggering transient interaction is found, then a possible next step is 
to see whether (in simulation) the presence of a PrPSc in the vicinity of a PrPC indeed would affect 
such an interaction. We take a two-step approach because a general exploration of possible transient 
interactions between PrPC and PrPSc by MD simulation would exceed our present computational 
capability, and because an accurate knowledge of the conformation of PrPSc is lacking. 

As mentioned in [1], a clear understanding of the physiological function of the PrPC and its 
interaction partner is still lacking, and this in turn prevents us from understanding the roles PrPC 
and PrPSc play in the pathogenesis of prion diseases. Thus, in the present work we institute 
structural changes in PrPC by hand, changes that we assume may be caused by hypothetical 
transient interactions, and follow the aftermath in each case by molecular dynamics (MD) simulation. 
Our goal is to investigate the dynamics of structural modifications of PrPC under simulations. 
Specifically we focus on the stability of HI after a change made in the S1-H1-S2 peptide (residues 
124-167). Our study contains two parts. In the first part we attempt to determine changes in which 
one of the two flanks - Fl and F2 - is more likely to initiate an unfolding of HI. We find the answer 
to be Fl. In the second part, we initiate specific structural changes in Fl, run MD simulations on the 
S1-H1-S2 peptide, and focus our attention on the way HI is affected. We find the native structure 
of the S1-H1-S2 peptide to be generally quite robust. Among structural alterations made to Fl, 
modification of the two dihedral angles of Asn^'*'^ is found to be most likely to lead to the unfolding 
of HI, and that when HI unfolds, it tends to form a /^-hairpin turn at residues 150-152, which is 
close to the 153-156 region reported in [50] ■ Our results also suggest that hydrophobic forces do not 
play a major role in the conversion process. 

2 Method 

2.1 Simulation parameter settings 

The simulation package AMBER 8 [6] is used for energy minimization and MD simulation with 
the AMBER force field ff03. As a prelude to each simulation a full conjugate gradient energy 
minimization is applied for 1000 iterations to allow the spatial positions of the atoms to relax to 
their respective local energy minima. During minimization and MD simulation SHAKE [32' is 
invoked to constrain hydrogen bonds. This has the effect of preventing the fast bond vibration 
motion of hydrogens. The cut-off distance for non-bonded interactions in all calculations is set to 
be 15 Angstrom. 

For the MD simulation, system temperature is at room temperature, or 300 K. Andersen tempera- 
ture coupling is employed to regulate temperature between protein and the environment. Simulation 
time step is 2 fs. Initial velocities of atoms in proteins are generated from Boltzmann distributions 
at 300 K temperature. The effect of solvent is represented by the modified generalized Born model 
of Onufriev et al. [28 , where the pH is set to neutral, and where, for calculating the effective Born 
radius, the maximum distance between a pair of atoms is set to be 12 Angstrom. Periodic boundary 
condition is not applied in the calculations. The cut-off distance of non-bonded interactions is set 
to be 15 Angstrom. 

The MD simulations are performed on a 32-node PC cluster at the National Center for High- 
performance Computing in Taiwan. Under the above settings the average CPU time needed to 
simulate the folding of a 40-amino-acid peptide for 1 ns is about half an hour. 

2.2 Designs for two series of simulations on S1-H1-S2 

Our intention is to represent the effect of potential conversion-triggering transient interactions on 
the S1-H1-S2 peptide (residues 124 to 167) by artificially induced structural changes to the peptide, 
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and study the stability of the affected peptide through MD simulation. Following this strategy we 
carry out two series of exploratory simulations classified according to the artificial changes made 
to the peptide and designated Enfm, where n= 1 and 2 is the classifier, and m enumerates the 
simulations in each class. In practice, a specific designation indicates a specific initial conformation 
for the peptide. In addition, fO (without the prefix En) denotes the benchmark simulation in which 
the S1-II1-S2 peptide has the native conformation as its initial conformation. 

The El series. The goal of this series of simulations is to identify which flank, Fl or F2, plays a 
crucial role in the stability of HI. We consider three extreme cases. In each case, either Fl, F2, or 
both together, is completely pruned from the S1-II1-S2 peptide and the remainder, initially in its 
respective native conformation, is simulated. In Elfl, F2 is pruned and the remainder consists of 
residues 124 to 154. In Elf2, Fl is pruned and the remainder consists of residues 142 to 167. In 
Elf3, both flanks are pruned and the remainder consists of residues 142 to 154. When Fl is pruned, 
native contacts of mouse PrPC including the hydrogen bonds in the anti-parallel /3-sheet, between 
the residue Arg^'^^ and the four residues Met^^*, Tyr^^^, Tyr^^^, and Pro^^*, and between residues 
Pro""^ and Tyr^^o are broken. In contrast no native contacts other than the hydrogen bonds in the 
anti-parallel /3-sheet are broken when F2 is pruned. It indicates Fl plays a more important role in 
the stability of HI than F2. This notion is discussed further in section [3l The total simulation time 
for this series is approximately 800 ns. 

The E2 series. Guilbert et al. [TS] pointed out that a major modification of the dihedral angles of a 
residue in Fl is required for the formation of a /3-sheet on the S1-H1-S2 peptide. This, together with 
the result from the El simulations, mainly that F2 plays a minor role in the stability of HI relative 
to Fl, motivate the E2 simulations describe below. In each simulation, the initial conformation of 
F2 relative to HI is unchanged, and one of two types of changes is made on Fl. In the first type, 
the dihedral angles of Asn^"*^, the residue joining Fl to HI, are changed. The native values of the 
dihedral angles are <I>o — —124.854° and '^q — 132.794°, which lie within the /3-strand region in the 
Ramachandran plot. In three simulations, designated E2fl, E2f2 and E2f3, the initial values of ($, 

are changed to ^')i = (60°, 60°), ^')2=($o, 0°), and ^')3=(-60°, ^o), respectively, see 
FigEl These modifications are made using DeepView/Swiss-Pdb viewer [M]. Figlll as is FigdJ are Fig. 2 is here, 
generated by Pymol f8|. In the Ramachandran plot, lies in the left-handed helical region, 

($, ^)2 in the right-handed helical region, and (<&, ^')3 in the /3-strand region. The total time for 
the three simulations is about 1 /iS. In all initial conformations in the above E2 series runs, native 
contacts including the hydrogen bonds in the anti-parallel /3-sheet, between the residue Arg^^^ and 
the four residues Met^^^, Tyr^^^, Tyr^^"^, Pro^^*, and between residues Pro^'^^ and Tyr^^° are broken. 
No new contacts are formed. 

In the second type, for which only one simulation (E2f4) is run, the initial internal conformation 
of Fl is altered by changing the values of dihedral angles of Ser^'^^ from $=-75,575° and ^'=150.389° 
to $=-75° and \l/=60°. This change has the effect of causing residues 124-143 to form a /3-hairpin-like 
structure (right panel in Figl2|). In this case, only the hydrogen bonds in the anti-parallel /3-sheet 
(in native mouse PrPC) are broken. 

2.3 In-house analysis tool 

We employ several in-house tools to analyze the simulation results. The code PTRAJ from the 
AMBER package is used to extract peptide conformation at 1 ns intervals, and the program DSSP 
[T9] to identify protein secondary structure. Two programs, g_sas and g_saltbr, from the MD pack- 
age GROMACS [38l [H |23] are used to calculate solvent accessible surfaces (SAS) and salt-bridge 
distances, respectively. We will use a similarity matrix to define a distance between two simulations. 
Confidence interval of the similarity thereafter is estimate using the maximum likelihood estimate. 
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3 Results 



3.1 Results of the fO simulation 

Stability of the S1-H1-S2 peptide. The benchmark experiment fO, was simulated for 200 ns 
during which, at intervals of 1 ns, the conformation of the S1-H1-S2 peptide is extracted. In the 
simulation of S1-H1-S2 (in mouse), the S1-H1-S2 forms a stable conformation different from its 
native state characterized by two quantities, the root mean square deviation (RMSD) and radius of 
gyration, shown in FiglH but HI remains largely intact. The RMSD calculates Ca atoms's positional 
difference between conformation in the simulation and the native conformation. 

Although the RMSD in the left panel of Fig. [3]shows two drops around 50 and 175 ns, the average 
RMSD is around 11 Angstrom during the whole simulation. The radius of gyration indicates spatial 
extent of all Cq atoms. The right panel of Fig. [3] shows that the spatial extent of the structure 
remains stable after collapsing around 12 ns. One may attribute the collapse to the following reason. 
Barducci et al.'s studies suggest Tyr^^^, Arg^^'', and Asp^^^ stabilize the /3-structure in the S1-H1-S2 
peptide [3]. Particularly, mutation of Asp^^* will severely cause disruption of the /3-sheet. Thus one 
should expect unzipping of the /3-sheet in our simulation because of absence of the Asp^^* in the 
initial conformation. As shown in Fig. [7] in next section there is the disruption of the /3-sheet after 
5 ns. In addition, there is a stable H-bonded turn around Met^'^* formed after 1 ns for the whole 
simulation. These two reasons may result in the collapse of conformation. Fig. 3 is here. 

Helical propensity in HI. The number of the residues from the HI region - D^'^'^WEDRYYR^^^ 
- forming the current a-helix is recorded at intervals of 1 ns in the 200 ns simulation. The result of 
the content variation along the simulation time is plotted in fO panel of Fig. 3) Since it takes 3.6 
residues to form a turn in an a-helix, one can define a zero-turn helix that consists of to 2 residues, 
one full turn helix contains residue 3 to 5, and two full turns helix contains 6 to 8 residues. This 
coarse-grained description provides a simple way to account for helical structural variations. In the 
fO simulation HI has zero turn in 8 accumulative ns out of a total 200 ns, one full turn in 70 out of 
200 ns, and two full turns in 122 out of 200 ns. There is a strong helical propensity in HI. These 
results are consistent with expectations 12^. Fig. 4 is here. 

3.2 Results of El simulations 

Stability of HI is more dependent on Fl. The number of residues forming the a-helix in HI 
region in Elfl, Elf2, and Elf3 is recorded. Each of the three simulations were run for 200ns. In each 
case, various data, including the number of residues forming the alpha-helix in HI region (residue 
144 to 152) were recollected at 50 ps intervals. The plots in Fig. [3]show the number of residues at 
every 20th data points, or effectively, data taken at 1 ns intervals. We will then apply a similarity 
analysis in the following to analyze the simulation results. 

First, similarity of the fO and the El simulations is defined as inner-product of two unit vectors, 
s. A vector characterized the simulation contains nine elements, which are accumulative times in 
200 ns simulation of 0, 1,..., 8-residue in HI forming a a-helix correspondingly. Then the unit 
vector is obtained by normalizing the vector by its length. The value of s ranges from to 1. A 
zero inner product indicates the two vectors to be completely different and one denotes the two 
to be identical. Thus one obtains a four by four similarity matrix to reveal all correlations among 
the fO and three El simulations S. Thereafter, we define a distance matrix, D = 1 — S to reveal 
similarities among the simulations. Fig. [5] shows an example of a schematic plot of family tree. 
Notes that each end point of branch in the tree represents a simulation and then total branch's 
length between two end points denotes their differences. For example length between the fO and the 
Elfl lfQ_Eifi = I fo-i+h-Eifi- Longer the branch's length is, the two simulations differ more. After 
twenty family trees corresponding to data recorded at different intervals are constructed, we apply 
the maximum likelihood estimate to evaluate confidence interval of length between two simulations. 
The 95 % confidence interval of length from the fO to the Elfl is Ifo^sifi = 0.1192 ± 0.0114 
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unit length. The lengths between the fO and two other simulations the Elf2 and the Elf3 are 
lfo-Eif2 = 0.172 ±0.0067 and Ifo-Eifi = 0.1875 ± 0.0102 respectively. 

The lengths calculations indicate the Ifo-Eifi to be the shortest at confidence level 0.95, which 
indicates Elfl to be most similar to fO, and the Elf2 and Elf3 are less similar to fO. Recall that in 
the Elfl simulation Fl is retained in the peptide, in the Elf2 simulation F2 is retained, and in the 
Elf3 simulation neither is retained, we thus conclude that relative to F2, Fl is significantly more 
crucial to the stability of the a-helix nature of HI and, by inference, that a transient interaction 
altering the structure of Fl is more likely to lead to a PrPC to PrPSc conversion than a transient 
interaction altering the structure of F2. Fig. 5 is here. 



3.3 Results of E2 simulations 

Native conformation of S1-H1-S2 has the lowest energy in simulation. Peptides in the 
E2 series have the same amino acid sequence but have different initial conformation. The total 
energy difference AE{fn) = EiE2{n)-E{m) = 41.7, 52.7, 37.9 and 23.2 kcal/mole for n= 1, 2, 3 and 
4 respectively, where for every case the energy is taken after the initial energy minimization and 
before the simulation begins. These results confirm the expectation that the native conformation of 
the S1-H1-S2 peptide has the lowest energy, at least compared to the initial conformations imposed 
on the peptide in the E2 series of simulations. This also provides a minimum necessary validation of 
the force field. We make a remark whose relevance will become clearer later: among the deformed 
peptides, E2f4 has the lowest initial energy. 

HI is unstable against modification in orientation of Fl. Recall that in simulations E2fl, 

E2f2 and E2f3, the initial relative orientation of Fl relative HI is changed (see Figl^]). FiglH] shows 

the number of residues that constitute the a-helix in HI as a function of simulation time in these 

simulations. Similarly, we apply the previous approaches to evaluate 95% confidence interval of Fig. 6 is here. 

length between the fO and four E2 simulations, which are lfa-E2fi = 0.1006 ± 0.0085, lfo-E2f2 = 

0.0933 ±0.0099, lfo-E2f3 = 0.0474 ±0.0062, and lfo-E2fi = 0.003 ±0.0014 respectively The results 

put the E2fl, E2f2, and E2f3 in an unstable class. The E2f4 and the fO are in the stable class. It 

indicates that retaining Fl in the peptide but changing its orientation relative to HI is sufficient to 

destabilize the structure of HI. In particular, the modification in E2fl, which has the longest length 

lfo-E2fi, has the largest effect on the destabilization of HI. 

HI is stable against modification in internal conformation of Fl. In simulation E2f4 the 
connection between Fl and HI is kept in its native state but the conformation of Fl is changed by 
altering the dihedral angles of Ser^^^ (see FigI5]). The simulation results are shown in the E2f4 panel 
of FiglH The shortest length between the fO and the E2f4, //0-B2/4 = 0.003 ± 0.0014, suggests that 
there is no difference between E2f4 and fO at 95% confidence level. The inference is that modifying 
the internal conformation of Fl does not destabilize HI. 

There are new hairpin-like turns in E2fl. Fig. [7] shows conformational transitions in fO, E2fl 

(representing E2fl, E2f2 and F2f3), and E2f4 during a 200 ns simulation, where "a-helix" includes 

the a-like structure aio-helix. In all three cases, a bending site giving a hairpin-like structure persists 

around residue 157 during the full simulation. Kozin et al. in their experimental studies on the sheep 

PrPC peptide in solution pointed out that this turn may form a part of the /3-sheet structure during 

conversion [20]. It is seen that fO not only retains helical structure most of time in the HI region Fig. 7 is here. 

(residues 144-152) but has a tendency to form additional helical structures around residues 126- 

134 as well. This last aspect is shared by E2f4. In contrast, E2fl rarely has any helical structure 

in the 126-134 region, has several more bends in the 132-144 region and, after 80 ns, many more 

hydrogen-bonded turns in 144-152 and a hairpin-like turn at residue 156. This hairpin-like turn is 

also identified to be part of the scrapie-like structure in the in silico study by Derreumaux ^llj . In 

addition, in the in vitro study by Megy et al., it was shown that the region 152-156, particularly at 

position 155, could be associated with the conversion [25], and that in PB buffer a mostly -hairpin 

like conformation was formed. Fig. 8 and ! 

are here. 
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The Glu^^^-Arg^'"'^ distance in E2fl is abnormally small. According to Kozin et al.'s studies, 
the overall stability of PrP142-166 is associated with the three salt bridges, bonds between two 
oppositely charged groups, in the HI region 20J. These bridges connect the pairs of residues Glu^"*^- 
^j.gi48^ ^gpi47_^j.gi5i^ g^j-^j Glu-'^^'^-Arg^^^. However, if one wants to stabilize an isolated HI helical 
structure, Dima and Thirumalai suggest the first and the third salt bridges should be Asp^'*^-Arg^^* 
and Asp^'**-Glu^^^ [12] for an isolated HI helix. Since our goal is to study overall stability of 
the peptide and only residues 145 and 167 of 1AG2 differ from the peptide they studied, we will 
still consider analyzing the three pairs suggested by Kozin et al.'s studies. The charges on the six 
residues involved in the ionic interactions are delicately balanced and the initial distances between 
two oppositely charged groups in the residue pairs for Glu^^^-Arg^''^, Asp^^^-Arg^^^, and Glu^^^- 
^j,gi56 jja^yg spans 1.052, 0.542, and 0.936 nm, respectively. These distances are measured from 
the NMR structure of the protein using the GROMACS tool g_saltbr, which computes the distance 
between the centers of mass of two oppositely charged groups. Notes that g_saltbr measures distance 
between center of mass of two oppositely charged groups. The left panel of FiglS] shows the distance 
of the Glu^"'^-Arg^'*^ as a function of simulation time in the three simulations fO, E2fl and E2f4. 
The histograms in the right panel gives the percentage of simulation time as a function of distance. 
The most likely distance in the fO and E2f4 simulations lies in the range 0.85 to 1.05 nm, whereas for 
E2fl the range is slightly greater, 0.95 to 1.1 nm. Thus the ionic bond in E2fl is slightly weakened 
in E2fl relative to the two other cases, while more rigidly confined in its range. Overall there is no 
significant difference among the three cases. The computed distances are broadly consistent with 
the NMR- measured distance of ~ 1.05 nm for the native conformation. 

Figini (right panel) shows the span of the Asp^^'^-Arg^^^ also do not vary much in the three 
simulations. In all case the span is concentrated within a relatively narrow range of 0.20 to 0.25 
nm. The computed value is however noticeably less than the NMR-measured distance of 0.54 nm 
for the native conformation. 

The situation shown in Fig. [10] for the Glu^^^-Arg^^^ is drastically different from the two other 
pairs. In the fO and E2f4 simulations the distance is mostly in the range 0.80 to 1.20 nm, and in the 
range 0.15 to 0.25 nm about 10% of the time. In sharp contrast, in the E2fl simulation, except for 
two transient periods, one at the beginning and one at around 70 ns, the distance is less than 0.25 
nm. That the distance settles to within a narrow range of 0.20 to 0.25 nm after 75 ns of simulation 
may be correlated with the (faint) appearance of a hairpin- like turn at residue 156 (see Fig. [7|). This 
distance is somewhat close to Megy et al. studies of sheep prion protein in PB buffer, which is about 
0.27 nm [25 . Although our simulations does not mimic PB buffer, it seems the modification of 
the relative orientation of the Fl flank to reach a hairpin-like conformation around residue 152-156 
similar to the Megy et al.'s discovery in PB buffer, which is believed to be likely associated with 
the pathogenic conversion. Namely, it may suggest the modification to be a required interaction for 
triggering conversation. It requires further investigations. Yet it is out of our current work's scope 
and will be left for future studies. Fig. 10 is here. 

Hydrophobic does not play a major role in the unfolding of HI. The solvent accessible 
surface (SAS) in a peptide gives an indication of effect of the hydrophobic force on its structure. In 
the mouse prion 1AG2 there are nine hydrophobic residues and five non-polar Glycines (underlined) 
in Fl: G^^^ LGGYMLG S AMSRPMIHFGN^^^: five hydrophobic residues in F2: N153mYRYPNQV 
YYRPVD167; and none in HI. 

The left panel in Fig[Tl] shows the SAS as a function of run time in the simulations fO, E2fl and 
E2f4. It is seen that in fO the SAS mostly fluctuates around 23 nni^ In E2fl, the SAS shows large 
fluctuation - between 18 and 26 nm^ - during the first 100 nm of the simulation, suggesting major 
conformational changes, but settles to a narrow range of 23±2 nm^ after 100 nm^. In E2f4, the 
fluctuation is larger then in fO, but does not have the large swings seen in E2fl. These fluctuations 
are only partly reflected in the histograms in the right panel, which gives percentage of time against 
SAS. The means and standard deviations in SAS for the three cases are 22.9±1.3, 23.0±2.0, and 
22.8±1.5 nm^ for fO, E2fl, and E2f4, respectively. The larger standard deviation in the SAS of 
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E2F1 is caused by the conformational fluctuation in E2fl during the early stages of the simulation. 
Overall, the SAS value does not appear to be sensitive to conformational transitions in the S1-H1-S2 
peptide. Fig. 11 is here. 



4 Summary and Discussion 

The template-assistance model attribute the pathogenesis of prion disease to an autocatalytic pro- 
cess, which occurs via transient interactions between PrPC (the mouse 1AG2 peptide) with PrPSc. 
Motivated by a basic assumption of the model, that a transient interaction may trigger a conforma- 
tion conversion, we carried out MD simulations of a simplified variant of the model. Specifically, we 
systematically imposed classes of simple alterations to the PrPC peptide, which served as represen- 
tations of the consequences of possible transient PrPC-PrPSc interactions, and examined whether 
such alterations would trigger the unfolding of the HI region of PrPC. 

The reliability of the force field used in the simulations (ff03 of the AMBER 8 package [5]) was 
verified in several ways: relative to its conformational variants, the S1-H1-S2 peptide had the lowest 
simulation energy in its native conformation; if the initial conformation of the peptide was native, 
then it would largely retain its native conformation during simulation; the computed spans of the 
three charged pairs of residues were in qualitative agreement with the measure values. The only 
quantitative exception was with the span of the Asp^'^^-Arg^^-'^; a computed distance of 0.20 to 0.25 
nm versus a NMR measured value of 0.54 nm. 

According to various in silico and in vitro studies [13l [TOl [37l [TlJ [25], the S1-H1-S2 peptide 
alone seems to be associated with the pathogenesis of the prion disease, and the disconnection of 
the H2-H3 segment from S1-H1-S2 seems to be a necessary condition for the latter's conversion to 
pathogenic form [TD]. All subsequent simulations were done on the S1-H1-S2 peptide composed of 
residues 124 to 163, the residues beyond F2 including the H2-H3 region - residues 164-226 - were 
left out. 

The first series of simulations on the S1-H1-S2 peptide were designed to show, of the two flanks 
Fl and F2, which would play a more important role in preserving the a-helical structure of HI. Four 
simulations were conducted: fO, the benchmark simulation on the entire peptide; Elfl, simulation 
on the peptide minus F2; Elf2, minus Fl; Elf3, minus both the Fl and F2. The similarity analysis 
of these results suggested that the integrity of the helical nature of HI depends crucially on the 
presence of Fl but only weakly on the presence of F2 (Fig. |4]). 

In the second series of tests F2 was left alone and artificial conformation alterations were made 
on Fl prior to simulation (Figl2|). In E2fl, E2f2, and E2f3, the relative orientations of HI and the 
Fl was modified by changing the two dihedral angles of residue Asn^''^ joining Fl to HI. In E2f4 
the relative orientations of HI and the Fl were unchanged while the internal native conformation 
of Fl was altered by changing the dihedral angles of Ser^'^^. It was found in the similarity analysis 
that keeping the Fl-Hl angle intact was crucial to the integrity of the native HI structure whereas 
keeping the internal structure of Fl was not (Fig. [6]). We remark that the present investigation is 
about the transient instability of HI after a disturbance is made to S1-H1-S2. Hence the simulation 
that follows is not supposed to be very long; in the present study the duration of all simulations were 
200 ns. Theoretically, a sufficiently long simulation will always bring the S1-H1-S2 peptide back to 
its native conformation, regardless of its initial state. 

Further examination of three other properties of the peptides - the existence of hairpin-like turns, 
the spans of charged pairs of residues, and the hydrophobic solvent accessible surface - showed results 
consistent with our interpretation that a (specific) change to the relative Fl-Hl orientation (E2fl, 
or E2f3, or to a lesser extent E2f2) would cause HI to unravel while a change in the internal 
conformation of Fl (E2f4) would not. The spans of all three pairs in the S1-H1-S2 peptide are 
similar during the fO and E2f4 simulations. In the E2fl, the span of the Glu^^^-Arg^^^ is reduces 
drastically from a native value of about 1.0 nm to about 0.2 nm. This change appears to correlate 



8 



with the appearance of an additional hairpin-hke turn around residue 152 in the simulation of E2fl 
(after 75 ns), a feature absent in the simulations of fO and E2f4. The SAS in the three simulations all 
average to about 23 nm^, but in the early part of the E2fl simulation (up to 100 ns) large fluctuation 
were seen, indicative of substantial conformational changes. 

The conformations of 1AG2 peptide and its subunits, the S1-H1-S2 peptide and the H2-H3 
domain, all turned out to be quite robust. Our simulations showed that neither the conformation 
integrity of S1-H1-S2 nor that of H2-H3 depends on the presence of the other. Furthermore, the 
native HI conformation was robust against any changes involving F2 and changes to the internal 
structure of Fl. This robustness is consistent with the fact that prion related diseases are not 
easily transmitted and rarely occur spontaneously, that is, non-infectiously. This may explain why, 
by and large, it afflicts only older people. Nevertheless, there does seem to be at least one type of 
vulnerability to this robustness: the helical structure of HI is prone to unraveling when the S1-H1-S2 
peptide suffers a large change in the relative Fl-Hl orientation. 
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Figure legends 



Figure [TJ Left: NMR structure of C-terminal of mouse PrPC; the peptides contains tlrree a-helices 
(H), two /3-strands (S), and six surface loops (L). Riglit: Tlic motif sequence, where the number of 
the leading residue in each motif is given. 

Figure [2] Initial conformations of the S1-H1-S2 peptide in the E2 series of simulations. Left: 
the initial conformations of E2fl, E2f2, and E2f3, where the native values of the dihedral angles 
Asn^'*^'s are changed. Right: the initial conformations of E2f4, where the dihedral angles of Ser^^^ 
are changed. 

Figure [3j Left panel plots root mean square deviation (RMSD) versus simulation time. Right panel 
plots radius of gyration. 

Figure |4| Stability of HI in El simulation. The number of residues in the a-helix as a function of 
simulation time in three El simulations. Result of fO is also shown for comparison. 
Figure \5\ A schematic plot of family tree for the fO and the El simulations. The branch length is 
scale free. 

Figure [6l Stability of HI in E2 simulations. Plots show number of residues in the a-helix as a 
function of simulation time in E2fl, E2f2, E2f3, and E2f4. 

Figure \7\ Conformational transitions of the 124-167 peptide in 200 ns simulations, with initial 
conformations being native (fO), E2fl, and E2f4, respectively. Y-axis denotes residue number and 
x-axis gives simulation time. Secondary structure are color-coded as shown at the top of the figure, 
where 3io-helix is classified as a-helix. 

Figure [8] Variation with time in the span of the Glu^'^^-Arg^'*^ in the fO, E2fl, and E2f4 simulations 

(left panel) and percentage of time the span has a specific distance (right panel). 

Figure [9] Variation with time in the span of the Asp"'^-Argi5i in the fO, E2fl, and E2f4 simulations 

(left panel) and percentage of time the span has a specific distance (right panel). 

Figure [TOl Variation with time in the span of the Glu^^^-Arg^^^ in the fO, E2fl, and E2f4 simulations 

(left panel) and percentage of time the span has a specific distance (right panel) . 

Figure [TT] Variation with time in the hydrophobic SAS in the fO, E2fl, and E2f4 simulations (left 

panel) and percentage of time the SAS has a specific value (right panel). 
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